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We analyze the transverse momentum spectra of charged pions measured in Au+Au 
collisions at = 200 GeV and in Pb+Pb collisions at ^/snn = 2.76 TeV using the 

Tsallis distribution modified to include transverse flow. All the spectra are well described 
by the modified Tsallis distribution in an extended transverse momentum range upto 6 
GeV/c. The kinetic freeze-out temperature (T), average transverse flow (/3) and degree 
of non thermalization (g) are obtained as a function of system size for both the energies. 
With increasing system size jd shows increasing trend whereas T remains constant. While 
the systems at RHIC and LHC energies show similar /3 and g, the parameter T is higher 
at LHC as compared to RHIC. The kinetic freeze-out temperature is also extracted using 
the measured charged particle multiplicity and HBT volume of the system as a function 
of system size and collision energies. 
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1. Introduction 

The Quantum Chromodynamics (QCD), the theory of strong interaction suggests 
that at energy density above ^ 1 GeV /ivp? the hadronic matter undergoes a phase 
transition to Quark Gluon Plasma (QGP)Pa phase where the relevant degrees of 
freedom are quarks and gluons. The heavy ion collisions at relativistic energy are 
the means to produce a large volume of hot/dense matter required to create and 
characterize such a phase.^^ 

The quark gluon matter presumably with local thermal equilibrium expands hy- 
drodynamically and undergoes a phase transition to hadronic matter which further 
cools till the multiple scatterings among particles are sufficient to keep them as one 
system. The hadrons then decouple from the system and their spectra would reflect 
the condition of the system at the time of freeze-out. Hadrons (pions, kaons and 
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protons) form the bulk of particles produced and are usually the first and easiest to 
be measured in a heavy ion collision experiment. Traditionally, statistical modeP 
has been used at SPS and RHIC energies to infer the conditions at freeze-out us¬ 
ing measured hadron ratios as input. Alternatively one can consider full transverse 
momentum (pr) spectra of hadrons in heavy ion collisions. The bulk and collective 
effect^^ show up in the low and intermediate pt regions of hadron spectra while 
the high pT region above 5 GeV/c consists of particles from jets which are produced 
in hard interactions. 

The Tsallis distributiorP® describes a system in terms of two parameters; tem¬ 
perature and q which measures deviation from thermal distribution. It has been 
shown in Refs.^^ that the functional form of the Tsallis distribution in terms of 
parameter q is the same as the QCD-inspired Hagedorn formulcf^^ in terms of 
power n. Both n and q are related and describe the power law tail of the hadron 
spectra coming from QCD hard scatterings. The Tsallis distribution has been used 
extensively to describe the pr spectra of identified charged hadrons measured in 
p + p collisions at RHIC and at LHC energies.^^^ It does not always provide the 
best description of hadron px distributions in heavy ion collisions which are mod- 
ihed due to collective flow and thus Tsallis blast wave method is used as in Ref.^ 
The average transverse flow can be included in Tsallis distribution and keeping the 
functional form to be analytical as done in Refs.^^^ The function presented in 
Ref.^ can be used in a wider px range as was done for both meson and baryon 
spectra for Au-|-Au collisions at ^snn = 200 GeV. 

We analyze the transverse momentum spectra of charged pions measured in 
heavy ion collisions. Recent measurements of identified charged particle spectra by 
PHENIX in different centralities of Au-|-Au collision at = 200 GeA^ and 

by ALICE in the most central (0-5%) and the most peripheral (60-80%) Pb-|-Pb 
collisions at ,/snn = 2.76 Te\^have been used in the study. The kinetic freeze-out 
temperature (T), average transverse flow (/3) and degree of non thermalization (q) 
are obtained as a function of system size for both the energies. As an alternative 
the (kinetic) freeze-out temperature is also extracted using the measured charged 
particle multiplicity and HBT volume of the system. 


2. Analysis of hadron spectra 

The transverse momentum spectra of hadrons can be described using the modified 
Tsallis distribution including the transverse flow as proposed in Ref.^ is given by : 


E 


(fN 

dp^ 


= Cn 




) —n 


( 1 ) 


Here is the normalization constant, nix = ^/Px i 7 “ — /32, (3 is the 

average transverse velocity of the system and T is the temperature. The power n 
is related to the non-extensivity parameter q as n = l/{q — 1). The parameter q 
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gives temperature fluctuationJi^ in the system as: 9 — 1 = Var{T)/{T)‘^. It can 
take a value between 1 and 4/3. Larger values of q correspond to smaller values of 
n which imply dominant hard QCD point-like scattering. Both n and q have been 
interchangeably used in Tsallis distribution.l^®2SH131 jn heavy ion collisions, the high 
PT tail decides the value of n. Phenomenological studies suggest that, for quark- 
quark point scattering, n ^ 4,1^3^ and when multiple scattering centers are involved 
n grows larger. When /3 is zero, Eq. [T]is the usual Tsallis equation which has been 
the most popular tool to characterize hadronic collision^i^^^^ in recent years. 
At low ptj Eq. [ 1 ] represents a thermalized system with collective flow and at high 
Pt it becomes a power law as follows 


E- 


dp^ 


~ Cn exp 


-j{mT - Ppt) 
T 


for Pt —t 0 , 




( 2 ) 


In this work, we focus on the study of the charged pion spectra measured in 
heavy-ion collisions at RHIC and LHC energies. The errors on the data are taken as 
quadratic sums of statistical and uncorrelated systematic errors. The RHIC mea¬ 
surements are available in pT range 0.5 - 6.0 GeV/c and we use the LHC measure¬ 
ments in the same range. The spectra are fitted with Eq. [T] and all the parameters 
are obtained as a function of system size (centrality) for both the energies. 

The freeze-out temperature T can also be extracted from the measured multi¬ 
plicity using following procedure. The particle number density n can be related to 
the measured particle multiplicity and HBT volume V as 


dN 

nV = —, 
dr] 


( 3 ) 


where dN/dr] is 1.5 times the total measured charged particle multiplicty {dNch/drf). 
The number density can also be expressed in terms of freeze-out temperature T 


n = Lfan{T)T'^. 


( 4 ) 


The parameter a„ (T) = gini{mi/T) where pi is the degeneracy factor and n- 
for meson species is given by 

1 n dx 


ii{mi/T) = 


2 X 1.2 Jq ^{y/x^ + (milT)'^ _ ^ 


( 5 ) 


The parameter a™ (T) = 3 for massless pion gas. In our study we assume that the 
system at freeze-out consists of pion {g = 3), kaon {g = A), p (g = 9), (p {g = 1), g 
[g = 1), ui {g = 3) mesons and obtained n/T^ as a function of temperature which 
is shown in Fig. [T] along with that for massless pion gas. 

The freeze-out temperature T can be obtained by numerically solving the fol¬ 
lowing equation 

1 1 dN 

(1.2/7r2)a„(r) 1/ d77 
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Fig. 1. (Color online) The variation of njT^ as a function of temperature for both hot meson gas 
(solid line) and massless pion gas (dashed line). 


The HBT volume V = R%de Riong, where Rside and Riong are the measured 

HBT radiiPl^ 

If transverse flow is present in the system, then the system volume obtained 
from measured HBT radii will be smaller than the fireball volume. To correct for 
this effect the HBT radii as a function of rriT are extrapolated to rriT = 0.140 
GeV/c^ which corresponds to px = 0. 


3. Results and Discussions 

Figure [5] shows the charged pion invariant yield spectra in Au+Au collisions at 
y/SNN = 200 GeV^ as a function of (mr - m), fitted with Eq. [T] The spectra 
in Au+Au collisions are given for 0-10%, 10-20%, 20-40%, 40-60% and 60-80% 
centralities, which are scaled up by factors given in the Fig. [2j Figure |3] shows the 
same for Pb+Pb collisions at -^/snn = 2.7(P^TeV fitted with Eq. [TJ It is seen that 
Eq. [1] describes the data in full px range for all collision centralities, both at RHIC 
and LHC energies. The collision centralities can be converted to average number 
of participants (Npart) using Glauber model which is proportional to initial system 
size. 

The parameter g as a function of (Npart) is shown in Fig. |4] (a). The value 
of q is higher for the peripheral collisions in comparision to other centralities, for 
both RHIC and LHC energies. The corresponding value for power n is shown in 
Fig.|4](b). The values of q (or n) are similar for RHIC and LHC which show similar 
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Fig. 2. (Color online) Invariant yield for charged pion^^as a function of {tut - ra) measured in 
Au+Au collisions at = 200 GeV for 0-10%, 10-20%, 20-40%, 40-60% and 60-92% centrality 

bins. The spectra are scaled up by a factor of 10, 5, 3, 2 and 1 for the respective centrality bins. 
The fitted Modified Tsallis function is shown by the black curve. 


degrees of thermalization at the two energy regimes. 

Figure[5](a) and[5](b) respectively show the kinetic freeze-out temperature (T) 
and the average transverse flow (/3) as a function of {Npart)- At RHIC energy, except 
for peripheral bin, the value of temperature remains constant for all centrality bins 
within uncertainties. For the most central collisions (0-5%), T has a higher value 
at LHC energy than that at RHIC. The transverse flow velocity /3 increases with 
the system size for both RHIC and LHC collisions and has almost same behavior 
in the two energy regimes. 

The freeze-out temperature is obtained from Eq. [6] using the measured particle 
multiplicity and HBT volume. Figure [6] (a) shows the freeze-out temperature as 
a function of system size and Fig. |6] (b) shows the same as a function of collision 
energy. The solid squares show the result obtained from the measured HBT radii. 
The open circles show the result obtained using the corrected HBT radii (at pr = 
0) as explained before. This correction makes the values of T slightly smaller. It 
is seen that the freeze-out temperature increases while going from RHIC energy to 
LHC energy. 

A comparison of Figs. [5] (a) and |6] (a) shows that the freeze-out temperature 
obtained from the two different methods follow similar trend as a function of sys¬ 
tem size. The temperature shows almost flat behavior with system size. Figures [5] 
(a) and |6] (b) show that the freeze-out temperature is more for LHC energy. The 
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Fig. 3. (Color online) Invariant yield for charged pionfl^l as a function of (rriT - m) measured in 
Pb+Pb collisions at = 2.76 TeV for 0-5% and 60-80% centrality bins. The Modified Tsallis 

fitted function is shown by the black curve. 
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Fig. 4. (Color online) Parameters (a) q and (b) n obtained from Eq. charged pions measured 
in Au+Au (open squares) collisions at = 200 GeV and in Pb+Pb (filled circles) collisions 

at y/SNN = 2.76 TeV as a function of (Npart)- 


temperature obtained from two measurements namely the pt spectra and the HBT 
measurements have same dependence on system size and energy. There is upto 20% 
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Fig. 5. (Color online) Parameters (a) T and (b) /3 obtained from Eq. [T]for charged pions measured 
in Au+Au (open squares) collisions at = 200 GeV and Pb+Pb (filled circles) collisions at 

yfSMN = 2.76 TeV as a function of (Npart)- 
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Fig. 6. (Color online) Kinetic freeze-out temperature obtained using Eq. [6] (a) as a function of 
(Npart) and (b) as a function of ,/sjfff (0-5% centrality) for RHIC and LHC energies. 


difference in the magnitudes which might arise due to the experimental error and 
different px range of the measurements affected by transverse flow differently. 

4. Conclusion 

The transverse momentum spectra of charged pions measured in Au-|-Au collisions 
at y/SNN = 200 GeV and Pb-bPb collisions at ^/snn = 2.76 TeV are analysed 
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using the modified Tsallis distribution. All the spectra used in this analysis are well 
described by this distribution. The parameter q of the modified Tsallis distribution 
suggests similar thermalization characteristics for systems at RHIC and LHC en¬ 
ergies. The kinetic freeze-out temperature extracted from pion px spectra remains 
flat for all centrality bins except for the peripheral bin at RHIC energy. The freeze- 
out temperature is higher at LHC energy than that at RHIC energy for the most 
central collisions. The transverse flow velocity increases with system size for both 
of these energies. The kinetic freeze-out temperature is obtained as a function of 
system size and collision energy, from the measurement of HBT radii and particle 
multiplicity. The measured HBT radii are extrapolated to mx = 0.140 GeV/c^ to 
correct for the effect of transverse flow. The freeze-out temperature obtained from 
particle spectra and HBT measurements show similar trend as a function of sys¬ 
tem size as well as a function of collision energies. However, the kinetic freeze-out 
temperature obtained using HBT radii remains larger than that obtained from the 
particle spectra. 
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